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Abstract 

We develop a graph-Hilbert-space framework, inspired by non-commutative 
geometry, on (infinite) graphs and use it to study spectral properies of 
graph- Laplacians and so-called graph-Dirac-operators. Putting the various 
pieces together we define a spectral triplet sharing most (if not all, depend- 
ing on the particular graph model) of the properties of what Connes calls 
a spectral triple. With the help of this scheme we derive an explicit expres- 
sion for the Connes- distance function on general graphs and prove both a 
variety of apriori estimates for it and calculate it for certain examples of 
graphs. As a possibly interesting aside, we show that the natural setting 
of approaching such problems may be the framework of (non-)linear pro- 
gramming or optimization. We compare our results (arrived at within our 
particular framework) with the results of other authors and show that the 
seeming diff^erences depend on the use of different graph-geometries and / or 
Dirac operators. 



1 Introduction 



We recently embarked on a programme to reconstruct continuous physics and/or 
mathematics from an underlying more primordial and basically discrete theory 
living on the Planck-scale (cf. As sort of a "spin-off" various problems 

of a more mathematical and technical flavor emerged which have an interest of 
their own. Discrete differential geometric concepts were dealt with in the 
theory of random graphs was a central theme of topics of dimension theory 
and fractal geometry were addressed in 0] . 

If one wants to recover the usual (differential) operators of continuum physics 
and mathematics by some sort of limiting process from their discrete protoforms, 
living on a relatively disordered discrete background like, say, a network, one has, 
in a first step, to study their discrete counterparts. This will be our main theme in 
the following with particular emphasis on discrete Laplacians and Dirac- operators 
on general graphs. 

We note in passing that functional analysis on graphs is both of interest 
in pure and applied mathematics and also in various fields of (mathematical) 
physics. For one, discrete systems have an increasing interest of their own or serve 
as easier to analyse prototypes of their continuum counterparts. To mention a 
few fields of applications: graph theory in general, analysis on (discrete) man- 
ifolds, lattice or discretized versions of physical models in statistical mechanics 
and quantum field theory, non- commutative geometry, networks, fractal geometry 
etc. From the vast and widely scattered literature we mention (possibly) very 
few sources which were of relevance for our own motivation or we came across 
recently (after finishing a first draft): to ||12|)111 and paper appeared 



after finishing our first draft, some more literature like e.g. ||Tj] was pointed out 



to us by Mueller-Hoissen; the possible relevance of references |T5| to [|T7| were 



brought to our notice by some unknown referee, [|18[ was also discovered by us 



only recently. Last, but not least, there is the vast field of discretized quantum 



gravitiy (see e.g. [|T9[ or [^). All this shows that the sort of discrete functional 
analysis we are dealing with in the following, is presently a very active field with 
a lot of different applications. 

We use the graph-Hilbert-space machinery, developed in the first part of our 
paper to investigate the spectral properties of graph Laplacians and Dirac oper- 
ators. In a next step we study and test concepts and ideas, which arose in the 
framework of non- commutative geometry. As we (and others) showed in preced- 
ing papers, networks and graphs may (or even should) be understood as examples 
of non- commutative spaces. A currently interesting topic in this field is the in- 
vestigation of certain distance functional on "nasty" or non-standard spaces and 
their mathematical or physical "naturalness" . Graphs carry, on the one hand, a 
natural metric structure given by a distance function d{x, y), with x, y two nodes 
of the graph (see the following sections). This fact was already employed by us 
in e.g. H to develop dimensional concepts on graphs. Having Connes' concept 
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of distance in noncommutative geometry in mind (cf. chapt. VI of |^), it is 
a natural question to try to compute it in model systems, which means in our 
context: arbitrary graphs, and compare it with the already existing notion of 
graph distance mentioned above. (We note in passing that the calculation of the 
Connes distance for general graphs turns out to be surprisingly complex and leads 
to perhaps unexpected connections to fields of mathematics like e.g. (non- jlinear 
programming or optimization; see the last section). 

Therefore, as one of many possible applications we construct a protoform of 
what Connes calls a spectral triple, that is, a Hilbert space structure , a corre- 
sponding representation of a certain (function) algebra and a (in our framework) 
natural candidate for a so-called Dirac operator (not to be confused with the 
ordinary Dirac operator of the Dirac equation), which encodes certain properties 
of the graph geometry. This will be done in section 4. 

In the last section, which deals with the distance concept deriving from this 
spectral triplet (as we like to call it), we will give this notion a closer inspection 
as far as graphs and similar spaces are concerned. In this connection some recent 
work should be mentioned, in which Connes' distance function was analyzed 
in certain simple models like e.g. one-dimensional lattices (^^-fS^). These 
papers already show that it is a touchy business to isolate "the" appropriate 
Dirac operator (after all, different Dirac operators are expected to lead to different 
geometries!) and that it is perhaps worthwhile to scrutinize the whole topic in a 
more systematic way. We show in particular that one may choose different Dirac- 
operators on graphs (or rather, different types of graphs over the same node set) 
which may lead to different results for e.g. the corresponding Connes-distance. 

The problem of finding suitable metrics on "non-standard" spaces is a particu- 
larly interesting research topic of its own, presently pursued by quite a few people 
(see the beautiful paper by Rieffel, and the references mentioned therein). 
Another earlier source is e.g. |T^. As to this latter paper we would like to remark 
that, while much of the framework is different from ours, there is, on the other 
side, a small overlap as far as some technical notions and results are concerned 
(after an appropriate translation of the respective technical notions and defini- 
tions). To give an example: While the definition of Dirac operators is different, 
some of the operator norms and metrics, being calculated, turn out to be identical 
to ours. This suggests a more careful comparison of the underlying conceptual 
ideas which we plan to give elsewhere. We presently extend this investigation 
of metric structures to lump spaces and probabilistic metric spaces (see and 



[0) in the general context of quantum gravity (cf. also p|). 

Our own approach provides a systematic recipe to calculate the Connes dis- 
tance in the most general cases of graphs and exhibits its role as a non-trivial 
constraint on certain function classes on graphs. We prove various rigorous a 
priori estimates and show how the constraints have to be dealt with in several 
examples. 
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Remark: (For possible reasons of priority) we would like to mention that many 
of our results can already be found in a (preliminary form) in an earlier draft 
version (p6|). 

2 A brief Survey of Differential Calculus on 
Graphs 

The following is a brief survey of certain concepts and tools needed in the further 
analysis. While our framework may deviate at various places from the ordinary 
one, employed in e.g. algebraic graph theory, this is mainly done for reasons 
of greater mathematical flexibility and generality and, on the other side, possi- 
ble physical applications (a case in point being the analysis of non- commutative 
spaces). Some more motivations are provided in |]l| and [^). We begin with 
the introduction of some graph theoretical concepts. We would however like to 
mention, that it is not our intention to cover any appreciable amount of the 
close interrelationship between graph spectra and graph characteristics (as has 
e.g. been done in ; our main emphasis lies on providing various Hilbert-space- 
techniques) . 

Definition 2.1 (Simple Locally Finite (Un)directed Graph) 

1. We write the simple graph as G := {V,E) where V is the countable set of 
nodes {rii} (or vertices) and E the set of bonds (edges). The graph is called 
simple if there do not exist elementary loops and multiple edges, in other 
words: each existing bond connects two different nodes and there exists at 
most one bond between two nodes. (We could of course also discuss more 
general graphs). Furthermore, for simplicity, we assume the graph to be 
connected, i.e. two arbitrary nodes can be connected by a sequence of con- 
secutive bonds called an edge sequence or walk. A minimal edge sequence, 
that is one with each intermediate node occurring only once, is called a path 
(note that these definitions may change from author to author). 

2. We assume the graph to be locally finite, that is, each node is incident with 
only a finite number of bonds. Sometimes it is useful to make the stronger 
assumption that this vertex degree, Vi, (number of bonds being incident with 
rii), is globally bounded away from oo. 

3. One can give the edges both an orientation and a direction (these two, in 
our view, slightly different geometric concepts are frequently intermixed in 
the literature) . In our context we adopt the following convention: If two 
nodes ni,nk are connected by a bond, we interpret this as follows: There 
exists a directed bond, dik, pointing from rii to rik and a directed bond, dki. 
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pointing in the opposite direction. In an algebraic sense, which will become 
clear below (for more details see also JJ/j, we call their superposition 



bik '■— dik — dki — —bki (1) 

the corresponding oriented bond (for obvious reasons; the directions are 
fixed while the orientation can change its sign). In a sense the above reflects 
the equivalence of an undirected graph with a directed multi-graph having 
two directed bonds pointing in opposite directions for each undirected bond. 

This way of algebraic implementation of geometric structures allows us to treat 
in principle all kinds of graphs on essentially the same footing. That is, it also 
applies to, say, graphs with only one directed edge being existent between two 
nodes. On the other side, an orientation should exist also for undirected graphs. 
As an aside, we remark that, on the one side, our generators, b^k, correspond to the 
oriented pairs of nodes, {i, k) in e.g. on the other hand, our dik correspond 
to the oriented pairs, (i, fc), in |T^. One sees from this that the conventions are 



far from being unique and that a certain unification may be desirable. 

We now take the elementary building blocks {ra,} and {dik} as basis elements 
of a certain hierarchy of vector spaces over, say, C with scalar product 

{ni\nk) = 6ik {dik\dim) = Su ■ 5km (2) 



Definition 2.2 (Vertex-, Edge-Space) The vector spaces (or modules) Cq, 
CI (a for antisymmetric) and Ci consist of the finite sums 

f-=^fini g-=^gikdik with gik = -gki and g' ■■=^gikdik (3) 

fiiQik ranging over a certain given field like e.g. C or ring like e.g. Z in case of 
a module. Evidently we have C Ci. 

These spaces can be easily completed to Hilbert spaces by assuming 

^|/,P<oo ^|c/,fcp<oo (4) 

if one chooses e.g. the field C (see the next section). Furthermore, one can 
continue this row of vector spaces in ways which are common practice in, say, 
algebraic topology ( see |l|] sections 3.1 and 3.2). In this context they are frequently 
called chain complexes (see also 0]). Evidently the above vector spaces could 
as well be viewed as discrete function spaces over the node-, bond set with rij, dik 
now representing the elementary indicator functions. 

In the same spirit we can now introduce two linear maps between Cq, Ci called 
for obvious reasons boundary- and coboundary map. On the basis elements they 
act as follows: 
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Definition 2.3 ((Co)boundary Operator) 

6 : dik ^ rik hence kk ^ nk - rn (5) 

d: rii^ ^{dki - dik) = ^ hi (6) 

k k 

and linearly extended. That is, 6 maps the directed bonds dik onto the terminal 
node and bik onto its (oriented) boundary, while d maps the node rii onto the sum 
of the ingoing directed bonds minus the sum of the outgoing directed bonds or on 
the sum o/ oriented ingoing bonds bki- 

The following results show, that these definitions lead in fact to a kind of 
discrete differential calculus on Cq,Ci. 

Observation 2.4 (Discrete Differential Forms) From the above it follows that 

df = d{^ fiUi) = ^(/fc - fi)dik (7) 

k,i 

Combining now the operators 5 and d, we can construct, what is called the 
canonical graph Laplacian. On the vertex space it reads: 

Observation 2.5 (Graph Laplacian) 

6df = - f^-"^^- = - - /*))^^ =■ -^f 

i k i k 

where Vi denotes the node degree or valency defined above and the k-sum extends 
over the nodes adjacent to rii . 



Note that there exist several variants in the literature (see e.g. |T5| or P|). Fur- 



thermore, many mathematicians employ a different sign-convention. We stick in 
the following to the convention being in use in the mathematical-physcis literature 
where —A is the positive(!) operator. 

This graph Laplacian is intimately connected with yet another important ob- 
ject, employed by graph-theorists, i.e. the adjacency matrix of a graph. 

Definition 2.6 (Adjacency Matrix) The entries aik of the adjacency matrix 
A have the value one if the nodes Ui^rtk are connected by a bond and are zero 
elsewhere. If the graph is undirected (but orientable; the case we mainly discuss), 
the relation between rii^rik is symmetric, i.e. 

aik = l aki = 1 etc. (9) 

This has the obvious consequence that in case the graph is simple and undirected, 
A is a symmetric matrix with zero diagonal elements. 
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Remark: More general A's occur if more general graphs are admitted (e.g. general 
multigraphs) . 

Observation 2.7 .■ With our definition of A it holds: 

A = A~V (10) 
where V is the diagonal degree matrix, having Vi as diagonal entries. 
(Note that the other sign-convention would lead to A = V — A). 

Proof: As we have not yet introduced the full Hilbert space machinery (which we 
will introduce in the next section), the proof has to be understood, for the time 
being, in an algebraic way. We then have: 

Af = A{J2f^n^) = J2MY,nk) = EE/'^^^ (^1) 

i k—i i k—i 

Vf = J2M)ni (12) 

i 

hence the result. □ 

(Here and in the following we use the abbreviation k — i if the nodes nk,ni are 
connected by a bond, the summation always extending over the first variable). 

As we already remarked above, our approach to functional analysis on graphs 
is perhaps a little bit different compared with the usual one. It seems therefore to 
be appropriate to exhibit some of the conceptual differences as compared to the 
more traditional framework, as e.g. developed in the beautiful monographs 
or 0, by briefly discussing the following (however only minor) point. In general, 
the graphs under discussion may be directed or undirected. In the traditional 
approach the edges are typically independently labelled of the nodes and the 
corresponding edge space, denoted in this case for the time being by Ci, is built 
over this edge set. In contrast to that habit we found it useful to label the 
occurring edges as di^, bik with bik = —bki which leads in our view to a more 
flexible discrete calculus and, among other things, to a natural Dirac operator on 
graphs (see below). 

More or less related to our operators d, d* {d* the adjoint of d; see the next 
section) are now the so-called incidence matrix, B, and its adjoint in the tradi- 
tional approach which relate the edges with the nodes. To do this, the edges are 
given an adhoc orientation, denoting one vertex arbitrarily as initial point, the 
other as end point. With the n labelled vertices, Ui and m labelled edges, Cj (Bij) 
has the entries 

+ 1 if nj is the endpoint of Cj 

— 1 if Ui is the initial point of Cj (13) 
otherwise 
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Evidently S is a mapping from C\ to Co and maps an edge to the respective 
difference of end vertex and initial vertex. By the same token, the transpose, 
is defined as a map from Co to Ci and one gets: 

B-B* -A (14) 

Note that the above introduced adhoc orientation does not enter in any end result; 
on the other hand, our approach is not based on such a contingent structure. 

3 Some Spectral Analysis and Operator Theory 
on (Infinite) Graphs 

After these preliminary remarks we now enter the heart of the matter. Our first 
task consists of endowing a general graph with a natural Hilbert space structure 
on which the various operators constructed in the following can operate. (The 
following analysis is done on undirected graphs, but could be extended to more 
general but less symmetric situations). 

Definition 3.1 (Hilbert Space) As indicated in the previous section, we ex- 
tend Co,Ci, Ci to the respective Hilbert spaces Hq, C Hi of sequences over 
the respective ON-bases {ui}, {dik}, that is {ni\nk) = Si^k, (c^ifeMi'fc') = ^n'^fefc'- 
As iJ", H we take the direct sums: 

H"" := Ho ®H^ CH := Ho © Hi (15) 

Note that members of if" can be written 

XI 9ikdik = 1/2 X Qikdik + 1/2 X gkidki = 1/2 ^ gik{dik - d^) = 1/2 ^ gikhk 

(16) 

Obviously H°- is a subspace of H and we have 

{hik\hik)^'2 (17) 

i.e. the bik are not(!) normahzed if the dik are. We could of course enforce this 
but then a factor two would enter elsewhere. 

With these definitions it is now possible to define the maps rf, 5 as true oper- 
ators between these Hilbert (sub)spaces. To avoid domain problems we assume 
from now on that the node degree v{ni) is uniformly bounded on the graph G, i.e. 

Vi < Vmax < OO (18) 
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Observation 3.2 We have the following relations 

d: Hq^H'^C. Hi, 6 : Hi ^ Hq (19) 

di^2 with 

di,2 ■ ^ 4i , ^ (20) 

respectively and linearly extended, are linear operators from Hq Hi and we 
have 

d = di-d2 (21) 

Similarly we may define 5 =: 5i and 62 via: 

(5i,2 : dik Uk, Hi (22) 

It is remarkable (but actually not surprising) that Vi < Vmax implies that 
all the above operators are hounded (in contrast to their continuous counterparts, 
which are typically unbounded). Taking this for granted at the moment, there are 
no domain problems and a straightforward analysis yields the following relations: 

Observation 3.3 

1. The adjoint d* of d with respect to the spaces Hq,Hi is 26 

2. On the other side we have for the natural extension of d, 6 to the larger 



space Hi (cf. the definitions in Observation \3.2J : 

Si = {diY , 62 = {d2y (23) 

hence 

(Si - ^2) = (di - d2y = d*^ 25 = 26i (24) 

3. Furthermore we have 

dl ■ di = 61 ■ di = d2 ■ d2 = V : Ui —>■ ViUi (25) 

d\ ■ d2 = 5i ■ d2 = d*2 ■ di = 82 ■ di = A : rii ^ Uk (26) 

k—i 

Similar geometric properties of the graph are encoded in the products coming 
in reversed order. 
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That and how d, d* encode some more geometric information about the graph 
can be seen from the following domain- and range-properties (for corresponding 
results in the more traditional approach see also p,p.24ff). 

Theorem 3.4 Let the graph he connected and finite, |V| = n, then 

dim{Rg{d*)) =n-l (27) 

dim{Ker{d*)) = ^Vi- (n-l) (28) 

i 

With dim{Hi) = ^jfi, dim{Hi) = 1/2 • dim{Hi) we have 

codim{Ker{d*)) = dim{Rg{d)) = n - 1 (29) 
We see that both Rg{d*) and Rg{d) have the same dimension {n — 1). 

Remark 3.5 In case the graph has, say, c components, the above results are 
altered in an obvious way; we have for example 

dim{Rg{d*)) = n - c (30) 

Proof: we first state the general result for bounded operators 

Rg{T*) = Ker{T)^ (31) 

we then have for T = d 

f e Ker{d) = d{f) = Y,ifk - mk (32) 

ik 

As the dik are linearly independent this entails fk = fi for the pairs (i.k) which 
occur in the sum. Since the graph is connected we have fk = fi = const for all 
nodes, hence dim{Ker{d)) = 1 and is spanned by X^i^j- ^^^^ proves the first 
item. 

In a similar way we proceed for d*. 

= d*{g) = "^C^igki - 9ik))ni ^(fi-fci - gik) = for all nodes (33) 

i k k 

In Hi gik,gki can be independently chosen. We have n linear equations, which 
are, however, not independent. There is, in fact, exactly one apriori constraint 
of the form 

Y.(T.(9^k - 9k^)) = (34) 

i k 

Hence, the above yields exactly n — 1 independent linear equations for the Yl '^i 
coefficients. This implies that the subspace, so defined, has dimension — 
{n — 1). This proves items two and three. □ 
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Observation 3.6 In the literature Ker{d*) is called (for obvious reasons) the 
cycle subspace (of e.g. j^). On the antisymmetric subspace we have d* = 26 
and 5{hik) = rik — rii. Choosing now a cycle, given by its sequence of consecutive 
vertices ra,^, . . . , rij^.; nj^.^^ := rijj, we have 

d*{J2^^.n^^) = 2 5^K+. -^0 = (35) 

that is, vectors of this kind lie in the kernel of d* 

We will now provide quantitative lower and upper bounds for the respective 
norms of the occurring operators. For d we have: 

i i k—i ik 

and it follows for the norm of the rhs: 

ik i k i^k 

= 2■Y,V^\f^?~2■Y,hf^ (37) 

i i^k 

The last expression can hence be written: 

IM/ir = m\Vf) - {f\Af)) = (/I - 2A/) (38) 

and shows the close relationship of the norm of d with the expectation values of 
the adjacency and degree matrix respectively the graph Laplacian. That is, norm 
estimates for, say, d, derive in a natural manner from the corresponding estimates 
for A or —A. It follows from the above that we have: 

Observation 3.7 

\\dfr = {f\d*df) = {f\-2Af) (39) 



I.e. 



and 



I.e. 



d*d = -2A hence \\df = sup (/| - 2A/) = || - 2A|| (40) 

11/11=1 



< sup (/I - 2A/) < 2vma. + 2 sup I < f\Af > I (41) 

11/11=1 11/11=1 



- A\\ < Vmax + \\A\\ (42) 
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We want to note that we are exclusively using the operator norm also for matrices 
(in contrast to most of the matrix literature), which is also called the spectral 
norm. It is unique in so far as it coincides with the so-called spectral radius (cf. 



e.g. or [HI), that is 

= sup{|A|; A e spectr{A)} (43) 



We now provide upper and lower bounds for the operator norm of the ad- 

we 



jacency matrix. A, both in the finite- and infinite-dimensional case. In p9 
estimated the upper bound by a method being different from the following calcu- 
lation which, being based on form- estimates, is much more direct. The previous 
proof was based on the so-called Gerschgorin-inequality for finite matrices and 
a not entirely straightforward extension to the infinite dimensional case. We ex- 
pect this upper bound to be well-known (cf. e.g. |T^, Theorem 2.8 - Lemma 
of Gabber- Gain - which is slightly more general). We do not know whether this 
is also the case for the lower bound. As such lower bounds are frequently less 
straightforward to derive and since, as a byproduct, we develop several potentially 
useful techniques, we give our own proof of the lower bound below. 

Theorem 3.8 (Norm of A) With the adjacency matrix A finite or infinite and 
a finite Vmax we have the following result ( a certain fixed labelling of the nodes 
being assumed): 

n 

limsup ■ ^^fj < ||^|| = sup{|A|; A G spectr{A)} < Vmax (44) 

i=l 

Proof: In order to prove the upper bound we use a form-estimate directly for the 
infinite case. We have 

|(x, Aa;)| = I xJa^^Xj \ < \xj\ ■ \xi\ (45) 

with a-'* = 1 or 0. Note that, due to the sjnnmetry of A, each term, \xj\ ■ \xi\ 
occurs twice in the above sum. With 

2\xj\ ■ \xi\ < + (46) 

we have 

rhs of (pD < t;^^,, ■ Ixp (47) 

and hence 

\{x,Ax)\ < v^axix,x) (48) 
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Strictly speaking we have the norm-bound up to now only proved for the above 
quadratic form. The well-known Riesz-lemma associates the form with a unique 
bounded operator which is the adjacency matrix we started with. This proves 
the first estimate. 

To prove the lower bound, we label the nodes or the corresponding orthonor- 
mal basis by (ei, 62, . . . ), and introduce the respective adjacency matrices on the 
corresponding subspaces, X„, belonging to the induced subgraphs, Gn, spanned 
by (ei, . . . ,e„). We choose a normalized vector, Xn, in X„ with all its entries 
being n~^/^. We then have 

n 

II All = sup I (X, Ax)\>\ AXn) I = I (Xn, A^Xn) I = ^l"^ ^ Vi (49) 

||x||=l ^ 

This proves the theorem. 

Lemma 3.9 The adjacency matrices, An, converge strongly to A and we have 
in particular \\An\\ y \\A\\. 



Remark 3.10 To prove strong convergence of operators is of some relevance 
for the limit behavior 0/ spectral properties of the operators An, A. That is (cf. 
e.g. section VIIL7), we have in that case (An, A selfadjoint and uniformly 
bounded) An —>■ A in strong resolvent sense, which implies that the spectrum of 
the limit operator. A, cannot suddenly expand, i.e. 

A G spec{A) ^ 3 A„ G spec{An) with A„ — A (50) 

and for a,b ^ speCpp{A) 

P{a,b){An) Pia,b)iA) strongly (51) 

Proof of the Lemma: Strong convergence can be proved as follows. A — An has 
the matrix representation 

^~ = ( j^t (J 1 (52) 



with 



where 



X 



Xr, 

Xrr- 



00 



with x„ = ^ XiCi , Xn = ^ XiCi (54) 



n+l 
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(and the Bn not to be confused with the incidence matrices of section 2) Multi- 
plying from the left with x we easily establish weak convergence since 

{x\{A - A 

and with n — cxd all the terms on the rhs go to zero, as for n — > oo 

since ||a;|| < oo and Bn,B^ are again uniformly bounded. 

To show strong convergence the critical term is B^^Xn- B^ maps the vector 
e Xn into — X Q Xn, Xn,X'^ living on the node sets Vn,V — Vn- As 

'^max < OO wc Can find for each given n a finite, minimal m„ so that all bonds 

beginning at nodes of end in V^^ , in other words: 

V n 3 m„ > n with e X^^ e X^ (56) 

or 

^m.n-'^rnn — ^^^{Xmn — Xn) (57) 

as Bl^^Xn = by construction. 

The Bi are uniformly bounded and ||2;to„ — — ^ foi" ^ — ^ oo, hence 

ll^m^a^mJI ^ with n ^ oo (58) 

Each I e N lies between some m„ and m(„+i) and we have 

B\xi = Bf{Xn + {Xl - Xn)) - Bf{xi - Xn) (59) 

as Bfxn — for all I > rrin > n. 

Z — > oo =^ m„ — > oo ^ n — > oo hence — — > (60) 
which shows that 

s - \im{Bfxi) = s - \im{Blx) = (61) 

□ 

Remark 3.11 A slightly simpler but perhaps less instructive proof can be given 
by exploiting the already established weak convergence together with special prop- 
erties ofAn,Bn etc., yielding 

\im{{A - An)x\{A - An)x) = - lim(A„x„|S„x;) = (62) 

n n 

To prove the monotone convergence of \\An\\ to \\A\\, we proceed as follows. 
For the principal minors we have 

An — PnAPn and An = PnAjnPn 

(63) 
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with Pn projecting on the subspace spanned by ei, . . . , e„ and m > n. Hence 

Pnll < and ||v4„|| < ||y4„|| (64) 

as \\Pn\\ = 1- From this we see that ||v4„|| is monotonely increasing with n ^ cxd 
and uniformly bounded by ||v4||. In other words: 

||A„||^a<p|| (65) 

The equahty of a and \\A\\ follows then immediately from the strong convergence 
of An towards A. This proves the above lemma. 

To test the effectiveness of the upper and lower bounds derived above, we 
apply them to a non-trivial model recently discussed in ||31|, i.e. the infinite binary 
tree with root Uq where Vq is two and Vi equals three for z 7^ 0. The authors show 
(among other things) that the spectrum consists of the interval [— 2a/2, 2v^], 
i.e. ||y4|| = 2^/2. Vmax IS three, we have to calculate limsupl/n ■ ^"fi. For 
simplicity we choose a subsequence so that n := n{N) with denoting the A^-th 
level (consisting of 2^ nodes) of the tree starting from the root hq. Note that in 
the corresponding induced subgraph the boundary nodes sitting in the A^-th 
level have only node degree one with respect to Gn but three viewed as nodes in 
the full tree. 

We then have 

JV n{N) N-1 N 

J]2\ ^1;, = 2 + 3-^2'= + 2^ = 3- ^2'= -2 -2^-1 (66) 

Hence 



n 

k=l i=0 k=l k=0 



n{N) N 

hm l/n{N) 5^ = 3 - 2 lim( J] 2^=-^)-^ = 2 (67) 
' i=o 



That is, our genral estimate imply 2 < \\A\\ < 3, which is not so bad. 



4 The Spectral Triplet on a general (undirected) 
Graph 

Note what we said at the beginning about our restriction to undirected graphs 
(made, however, only for convenience!). Furthermore our Dirac operator inter- 
twines node-vectors and bond-vectors while in other examples it maps node- 
to node-functions. Our bond-functions have (in some sense) the character of 
cotangential-vectors, while in other approaches derivatives of functions are inter- 
preted as tangent-vectors. In our view, the latter formalism is effective only in 
certain classes of highly regular models (like e.g. lattices) where one has kind of 
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global directions and will become cumbersome for general graphs. We developed 
this latter approach a little bit in section 3.3 of [|l| and showed how these cotan- 
gent and tangent vectors can be mapped into each other. 

The Hilbert space under discussion in the following is 

H = Ho®Hi (68) 
The natural representation of the function algebra 

{/;/GCo,sup|/,|<oo} (69) 

i 

on H by bounded operators is given by: 

^o: f-f = J2M-n,foTf'eHo (70) 



Hi : f ■ '^^dikdik '■= ''^fiQikdik (71) 



From previous work ([0]) we know that Ci carries also a right-module structure, 
given by: 

dikdik ■ f '■= Qikfk ■ dik (72) 

(For convenience we do not distinguish notationally between elements of and 
their Hilbert space representations). 

An important object in various areas of modern analysis on manifolds or in 
Connes' approach to noncommutative geometry is the so-called Dirac operator D 
(or rather, a certain version or variant of its classical counterpart; for the wider 
context see e.g. |^ or |^ to [^). As D we will take in our context the operator: 

--(It) (-) 



acting on 



Ho 



with 

d* = (5i - 62) (75) 

Note however, that there may exist in general several possibilities to choose such 
an operator. On the other hand, we consider our personal choice to be very 
natural from a geometrical point of view. 
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Lemma 4.1 There exists in our scheme a natural chirality- or grading operator, 
X and an antilinear involution, J. given by 



with 

[X,^] = x-D + D-x = (77) 



and 



so that 



J-f-J = f (79) 



These are some of the ingredients which estabhsh what Connes calls a spectral 
triple (cf. e.g. [0] or [^). We do not want, however, to introduce the full 



machinery at the moment as our scheme has an independent geometric meaning 
of its own. So, being careful, we call in the following these structures simply 
spectral triplets (we were kindly warned by B.Iochum to be more careful with this 
concept; note the observation below about the non-compactness of the inverse of 
such Dirac-operators on infinite graphs with a uniformly bounded vertex degree). 



Definition 4.2 (Spectral Triplets) ^45 spectral triplet on a general (undirected) 
graph we take 

{H,J^,D) (80) 

At this point we would like to remark the following. In our general framework 
we restricted ourselves, mostly for (possibly subjective) aesthetic reasons - the 
mathematics tends to be more transparent - to undirected graphs and a total 
Hilbert space being the direct sum of the node space (a function space) and 
the bond space (sort of cotangent vectors). A was then selfadjoint and a Dirac 
operator emerged naturally as kind of a square root of the Laplacian. 



On the other side, if one studies simple models as e.g. in p3| to 25 , other 



choices are possible. In where the one-dimensional lattice was studied. 



the symmetric difference operator was taken as Dirac operator. In the one- 
dimensional lattice was assumed to be directed (i.e. only were present) 
and the Dirac operator was defined as a certain self adjoint "doubling" of the 
(one-sided, i.e. non-symmetric) adjacency matrix. This latter model would fit 
in our general approach if we had included more general graphs. All these Dirac 
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operators are different and it is fience no wonder tfiat tfiey lead to different 
consequences (sec below). It is our opinion that, in the end, an appropriate 
choice has to be dictated by physical intuition. Nevertheless, this apparent non- 
uniqueness should be studied more carefully. 

As can be seen from the above, the connection with the graph Laplacian is 
relatively close since: 



and 



d*d = -2A (82) 

dd* is the corresponding object on Hi. (In the vector analysis of the continuum 
the two entries correspond to divgrad , graddiv respectively ) . 

Observation 4.3 Note that all our operators are hounded, the Hilbert space is 
(in general) infinite dimensional, hence there is no chance to have e.g. {D — z)~^ 
or {D^ — z)"^ compact. At the moment we are sceptical whether this latter phe- 
nomenon dissappears generically if the vertex degree is allowed to become infinite. 
There are some results on spectra of random graphs which seem to have a certain 
bearing on this problem. 

We now calculate the commutator [D, f] applied to an element f & Hq: 

{d-f)f^J2{fkfl-fJl)d,, (83) 

ik 

{f-d)f^Y.M-f'^d,, (84) 

ik 

hence 

[D,f]f^Y.^fk-fi)fkdik (85) 

ik 

On the other side the right-module structure allows us to define df as an operator 
on Hq via: 

df-f'= [J^Uk - fi)d^k) ■ (5^ flnu) = Y^Uk - fi)f'kd,k (86) 

ik k ik 

In a next step we define df as operator on Hi which is not as natural as on Hq. 
We define: 

df\H^ ■ dik {fi - fk)nk (87) 
and linearly extended. A short calculation shows 

df\H, = -{df\HS = [d\n (88) 
This then has the following desirable consequence: 
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Observation 4.4 With the above definitions the representation of df on H is 
given by 

and it immediately follows 

5 The Connes-Distance Function on Graphs 

From the general theory we know that: 

||r|| = ||r*|| (91) 

Hence 

Lemma 5.1 

\\[dJ]\\ = mf]\\ = \\[d*J]\\ (92) 

and 

ll[A/]|| = ||M,/]|| (93) 



Proof: The left part of (|92D is shown below and is a consequence of formula 



the right identity follows from (pi|). With 

X:-{;) (94) 

and Ti := [rf, /], T2 := [d*, /], the norm of [D, f] is: 

\\[D, f]f = sup{||Tixf + WT^yf; ||xf + \\yf = 1} (95) 

Normalizing now x,y to ||x|| = \\y\\ = 1 and representing a general normalized 
vector X as: 

X = \x + fxy , A, /i > and + /i^ = 1 (96) 

we get: 

II [D, f]f = sup{Alrix|r + AT2yf; \\x\\ = \\y\\ = 1, A^ + /i^ = 1} (97) 
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where now x, y can be varied independently of A, /i in their respective admissible 
sets, hence: 



\[D,f]f = snp{X^T,f + fi'\\T2f} = \\T,f □ (98) 



(as a consequence of equation (p2D). 

It follows that in calculating ||[Z}, /]|| one can restrict oneself to the easier to 
handle /]||. For the latter expression we then get from the above {x G Ho): 

Vj 

i k=l 

Abbreviating 

Vj 

J2\fk-fi\'=:a,>0 (100) 
k=l 

and calling the supremum over i a^, it follows: 

\\df ■ = as ■ C^ai/tts ■ |xip) < (101) 

i 

for = J2i kip = 1- 

On the other side, choosing an appropriate sequence of normalized basis vec- 
tors so that the corresponding converge to we get: 

IM/-e,||2^a, (102) 

We hence have 
Theorem 5.2 

\\[DJ]\\=snp{J2\fk-m'/' (103) 

' k=i 

The Cannes- distance functional between two nodes, n, n', is now defined as 
follows: 

Definition 5.3 (Connes-distance function) 

dzstcin,n') := sup{|/„, - \\[D,f]\\ = \\df\\ < 1} (104) 



Remark 5.4 It is easy to prove that this defines a metric on the graph. 
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Corollary 5.5 It is sufficient to vary only over the set {/; \\df\\ = 1}. 



Proof: This follows from 

\fk-fi\=c-\h/c-f,/c\; c=\\df\\ (105) 

and 

\\d{f/c)\\=c-^df\\=l (106) 
with c < 1 in our case. □ 

It turns out to be a nontrivial task (in general) to calculate this distance on 
an arbitrary graph as the above constraint is quite subtle . The underlying reason 
is that the constraint is, in some sense, inherently non-local. As / is a function, 
fn' — fn has to be the same independently of the path connecting n' and n. On the 
other side, in a typical optimization process one deals with the individual jumps, 
fk — fi, along some path. It is then not at all clear that these special choices can 
be extended to a global function without violating the overall constraint on the 



expression in theorem |5.2| . Nevertheless we think the above closed form is a solid 
starting point for the calculation of distc on various classes of graphs or lattices. 
We discuss two examples below but refrain at this place from a more complete 
treatment, adding only some observations concerning the relation to the ordinary 
(combinatorial) distance function introduced in the beginning of the paper. 

Having an admissible function / so that supj(^^!^]^ \ fk — fil'^Y^'^ < 1, this 
implies that, taking a minimal path 7 from, say, n to n', the jumps \fu+i — fu\ 
between neighboring nodes along the path have to fulfill: 

\fu+i-fu\<l (107) 

and are typically strictly smaller than 1 as long as there are not a sufficient 
number of "zero-jumps" ending at the same node. 

On the other side the Connes distance would only become identical to the 
ordinary distance d{n, n') if there exist a sequence of admissible node functions 
with all these jumps approaching the value 1 along such a path, which is however 
impossible in general as can be seen from the structure of the constraint on the 
expression in theorem ^]2| . Only in this case one may have a chance to get: 



I - /.) I ^ 5^ 1 = length{^) (108) 

7 7 

We express this observation in the following way 

Observation 5.6 (Connes-distance) One has within our general scheme the 
following inequality 

distc{n,n') < d{n,n') (109) 
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A fortiori one can prove that distc between two nodes in an arbitrary graph is even 
smaller than or equal to the corresponding Connes- distance taken with respect to 
the (one- dimensional) sub- graph formed by a minimal path between these nodes, 
i.e. 

distc{n,n') < distc{rnin.path){n,n') (HO) 

The simple reason is that one has more admissible functions at ones disposal 
for a subgraph, hence the supremum may become larger. This latter distance, 
on the other side, can be rigorously calculated (see Example 2 below) and is for 
non-neighboring nodes markedly smaller than the ordinary distance. 



Corollary 5.7 The last inequality implies also that with G' an induced subgraph 
ofG It holds (n,n' e V C V): 

distcin, n'] G') < distc{n, n'; G) (111) 

We remarked above that the calculation of the Connes distance on graphs 
is to a large part a continuation problem for admissible functions, defined on 
subgraphs. Then the following question poses itself. 

Problem 5.8 For what classes of graphs and/or subgraphs do we have an equality 
in the above corollary? 



Remark 5.9 Equality can e.g. be achieved for trees. Other results in this direc- 
tion are in preparation. 



These general results should be contrasted with the results in to [ p5[] . 
Choosing e.g. the symmetric difference operator as Dirac operator in the case of 
the one-dimensional lattice the authors got in ^ a distance which is strictly 
greater than the ordinary distance but their choice does not fulfill the above 



natural constraint given in Theorem 5.2. Note in particular that our operator 



d is a map from node- to bond-functions which is not the case in the other 



examples. In |]25[ the authors employed a symmetric doubling of the upper half 
of our symmetric adjacency matrix as Dirac operator. In the case of the one- 
dimensional (directed) lattice this then leads (so to say) to only one (directed) 
bond per node and makes the optimization process quite simple, hence leading 
to the ordinary distance which would have been also the case in our general 
scheme had we admitted directed graphs. We conjecture however that for more 
general graphs a relation related to the one given in Theorem ^. 2| would enforce 
the Connes-distance to be again strictly smaller than the ordinary distance for 
non-neighboring points. This is however an interesting point and we plan to 
discuss generalisations of our framework and more general examples elsewhere. 
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We want to close this paper with the discussion of two examples. The first 
one is a simple warm-up exercise, the second one is the one-dimensional lattice 
discussed also by the other authors mentioned above (treated however within their 
respective schemes) and is not so simple. The technique used in approaching the 
second problem may be interseting in general. While we solved it starting, so to 
speak, from first principles, the real mathematical context, to which the strategy 
is belonging, is the field of (non-)linear programming or optimization (see e.g. 
35] or any other related textbook). This can be inferred from the structure of 



our constraint on the expression in theorem ^.2| . This means that the techniques 



developed in this field may be of use in solving such quite intricate problems. 

Example 1: The square with vertices and edges: 

Xi — X2 — X3 — X4 — Xi (112) 

Let us calculate the Connes-distance between xi and 0:3. As the sup is taken 
over functions(!) the summation over elementary jumps is (or rather: has to be) 
pathindependent (this is in fact both a subtle and crucial constraint for practical 
calculations). It is an easy exercise to see that the sup can be found in the class 
where the two paths between xi, X3 have the valuations (1 > a > 0): 

Xi- X2: a , X2- X3 : (1 - a^y^'^ (113) 

Xi — X4 : (1 — a'^y^'^ , Xi — X3 : a (114) 

Hence one has to find supo<„<]^(a + Vl — ct^)- Setting the derivative with respect 
to a to zero one gets a = ^/l/2. Hence: 



Example 5.10 (Connes-distance on a square) 

distc{xi, X3) = V2 < 2 = d{xi, X3) (115) 
Example 2: The undirected one- dimensional lattice: 

The nodes are numbered by Z. We want to calculate distc{0,n) within our 
general framework. The calculation will be done in two main steps. In the first 
part we make the (in principle quite complicated) optimization process more 
accessible. For the sake of brevity we state without proof that it is sufficient to 
discuss real monotonely increasing functions with 

„ / , N I I \ - ; for k < , ^ 

fik) ={ r) [ r ~ (116) 

for k > n 
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and we write 

k 

f(k) = /(O) + 0<k<nhi>0 (117) 

i=l 

The above optimization process then reads: 

Observation 5.11 Fine? sup ^"^^^ /ij under the constraint 

hl<l,hl + hl<l,... ,hl + hl_, <l,hl<l (118) 

The simphfying idea is now the following. Let h :— (/ii)f=i be an admissible 
sequence with all h^_^i + < 1. We can then find another admissible sequence 
h' with 

J2K>J2h^ (119) 

Hence the supremum cannot be taken on the interior. We conclude that at least 
some hf_^i + /i? have to be one. There is then a minimal i for which this holds. 
We can convince ourselves that the process can now be repeated for the substring 
ending at i + 1. Repeating the argument we can fill up all the entries up to place 
i + 1 with the condition hf_^i + hf — 1 and proceeding now upwards we end up 
with 

Lemma 5.12 The above supremum is assumed within the subset 

hl<l,hl + hl = l,... , hl_, + hl = l,hl<l (120) 

This concludes the first step. 

In the second step we calculate sup |/(0) — f{n) \ on this restricted set. From 
the above we now have the constraint: 

hl<l,hl^l-hl,hl^hl,hl^l-hl,... ,hl^l- hi or hi (121) 

depending on n being even or uneven. This yields 

{1 forn = 1 

(n/2) • sup(/ii + v"! - h\) = (n/2) • ^2 for n even 
sup([n/2] • {hi + y^l — hi) + hi) for n uneven 

(122) 

In the even case the rhs can be written as \/n?/2 = ^[n^/2]. In the uneven case 
we get by differentiating the rhs and setting it to zero: 

/ir" = AJ,/l + Al , Vl - {hrn' = 1/^1+^ (123) 
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with An = 1 + l/[n/2]. We see that for increasing n both terms approach l/-\/2, 
the result in the even case. Furthermore we see that the distance is monotonely 
increasing with n as should be the case for a distance. This yields in the uneven 
case 

which is a little bit nasty. Both expressions can however be written in a more 
elegant and unified way (this was a conjecture by W.Kunhardt, inferred from 
numerical examples). For n uneven a short calculation yields 

[n^/2] = (n^ - l)/2 = 1/2 • (n - l)(n + 1) = 2[n/2]([n/2] + 1) (125) 

(with the floor- , ceiling-notsition the expressions would become even more ele- 
gant). With the help of the latter formula the rhs in ( 124 ) can be transformed 
into 

rhs of (|124D = yp/2[TT (126) 



Conclusion 5.13 For the one- dimensional undirected lattice we have 

I \/\n^ 1% for n even 
distci^^n) = / \ (127) 

I V F /^-l + ^ J'^''^ ^ uneven 

Remark: With the help of the methods, introduced above, we can now estimate 
or rigorously calculate the Connes-distance for other classes of graphs. 

Acknowledgement: We thank the referees for their constructive criticism. 
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